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Abstract 

We study interactions of electro-magnetic fields with the curvature 
tensor of the form \R^ va pF^ v 'F a @ . Such coupling terms though are 
invariant under general coordinate transformation and CPT, however 
violate the Einstein equivalence principle. These couplings do not 
cause any energy dependent dispersion of photons but they exhibit 
birefringence. We put constraints on the coupling constant A using 
results from solar system radar ranging experiments and millisecond- 
pulsar observations. We find that the most stringent constraint comes 
from pulsar observations and is given by A < 10 11 cm 2 obtained from 
the timing of binary pulsar PSR B1534+12.. 



1 Introduction 

Einstein's general relativity is based on the twin principles of general co- 
ordinate invariance and the Equivalence principle. A generalization of the 
standard minimal coupling of electro-magnetic fields in curved space-time 
introduced in [1] 

£ = ~\ F M „ + A R>»»» F^F pa (1) 



is a theory which violates the Einstein's equivalence principle (which 
states that in a local inertial frame the theory should reduce to the form 
given by special relativity) as the Riemann curvature tensor does not van- 
ish in a local inertial frame. This theory however is invariant under general 
coordinate transformation as well as CPT. 

In this paper we study some of the experimental consequences of the 
curvature couplings of electro-magnetic fields and put constraints on the 
coupling constant A from radar ranging past the sun [2] and from the ob- 
servations of millisecond pulsars [3]. The effect of the curvature term in the 
lagrangian is to make the photon orbits polarization dependent. Therefore 
the two independent polarizations of the electromagnic waves experience dif- 
ferent amounts of bending and time delay in a gravitational field. There is 
however no energy dependent dispersion of the photon trajectories. 

Constraints on the dispersion relations of the generic form p 2 — E 2 = 
(E/M) n can be obtained from the non observation of time lag between GRB 
signals of different energies which have propagated over cosmological dis- 
tances. Bounds of the order of M > 10 l5 gev [6] have been obtained from 
the non-observation of an energy dependent time delay of Gamma Ray Burst 
signals. 

An example of a coupling term which does not give rise to an energy 
dependent group velocity is given by Lj = A K pvpu F^ v F pa [7], where K pupa 
is a non-dynamical constant which has the same symmetries under exchange 
of indices as the Riemann tensor. This coupling violates CPT and Lorentz 
invariance and it gives rise to optical rotation of photon signals. Kostelecky 
and Mewes [7] put upper bounds ~ 3 x 10~ 32 on the components of K from 
the non-observation of optical activity in signals from radio galaxies. 

In the present example we find that due to the Riemann term, the cou- 
pling becomes vanishingly small for cosmological background metrics and the 
advantage of large cosmological distances over which the effect of dispersion 
or polarisation is accumulated in the examples discussed earlier are lost. We 
find that the strongest bounds come from the consideration of compact stars. 
Radar echo experiments of radio signals passing in the vicinity of the sun give 
a bound of A < 3.9 x 10 19 cm 2 . This is two orders of magnitude more stringent 
than bounds on A obtain from the birefringent bending of light by the sun 
[8]. We find that the most stringent bound on A is obtained from the timing 
of binary pulsar PSR B1534+12 signals [3] and is given by A < 0.6 x 10 n cra 2 . 
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2 Curvature coupling electrodynamics 

The equation of motion for F^ v from the Lagrangian (1) is given by 

V„F^ = 2A [ir^(V„F pCT ) + [VvR^Fpv] (2) 
In addition we have the Bianchi identities 

V^V + VaR^ffp + V P R»\ P = (3) 

Setting P — v in (3) we get 

V„iTV = V.i^ - V p i2" ff (4) 

Using (4) in the equation of motion (2) we obtain 

V V F» V = 2A {R^{V u F pa ) + (V P R£)F pa ) (5) 

In the absence of sources (R^ = 0) equation (5) simplifies to 

V„F^ = 2A R^{V v F pa ) (6) 

To obtain the trajectory of photons in curved space it is convenient to 
choose a locally flat inertial frame and derive the dispersion relations for the 
photon momenta in the inertial frame. The experimentally observed 
photon momenta in the coordinate frame k^ is then obtained by making use 
of the tetrads k^ = e^p^y In our application the photon wavelengths will be 
much smaller than the curvature scale of the gravitating bodies and therefore 
we can make the eikonal approximation 

V M FMM=p M FMM (7) 

in which case the equations of motion (6) in a local inertial frame appear as 

p {v) F<")<"> = 2\R^^p (l/) F {am (8) 

We also have the Bianchi identity, 

P{a)F(p)(y) + P(is)F(a)(p) + P(p)F (l/) ( a) = (9) 
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Setting fx — j in (8) ( we denote the time index by and spatial indices by 
Latin and four indices by Greek letters) we obtain 

p {0 )F^°) +p w F®W - 2A(^)MWWp M F (a)(0) + R mu){a){k) P(u)F {am ) = 0(10) 

We use the Bianchi identity (9) to write the magnetic field tensor in terms 
of electric field tensor 

P(o) F (j)(k) = P(k)F(j)(p) - P(j)F(k)(0) (11) 

Using (11) in (10) we finally obtain the wave equation for the electric field 
polarization vector in a local inertial frame, 

(p^poo^j - p (0) m + ^'^(^w) F(km = (12) 

Using (8) to substitute for the second term we find that (12) reduces to the 
form 

(?) 

(P^WgJ + 4A(-^ e (0) (fc) + = (13) 

P(o) 

where 

^w^'^Wjpm ( 14 ) 

The dispersion relations for electromagnetic waves can be obtained by setting 

Det\pV>p w 6$ + 4A(-^e% + e%)] = (15) 

P(o) 

For the Schwarzschild metric 

rfs 2 = (1 _ 2^ )rft 2 _ (1 _ 2 ^)-l rfr 2 _ r 2 {dQ 2 + Sin 2 dd(j) 2^ 

r r 
the non-zero components of Riemann tensor in a local inertial frame are 

2GM 

«(0)(1)(0)(1) = -«(2)(3)(2)(3) = " 

fl(0)(2)(0)(2) = #(0)(3)(0)(3) = --R(l)(2)(l)(2) = --R(l)(3)(l)(3) = (17) 
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and the components of e( a )(/3) are 



C(0)(0) = 


A/T 

r3 ( 2 P(1) 2 +p (2) 2 +p (3) 2 ) 


e (l)(l) = 


r 3 ( 2 P(0) +P(2) +P(3) ) 


6(2)(2) = 


r^A/r 

r 3 (P(0) P(l) + 2p (3 ) ) 


C(3)(3) = 


KjrlVl Q 2 n 2\ 
r 3 (P(0) P(l) + 2p (2 ) ) 


C(0)(1) = 


2GM 

r 3 P(0)P(1) 


C(0)(2) = 


CM 

r 3 P(0)P(2) 


C(0)(3) = 


GM 

r 3 P(0)P(3) 


£(1)(2) = 


GM 

r 3 P(1)P(2) 


C (l)(3) = 


r 3 P(1)P(3) 


C(2)(3) = 


2GM 

r 3 P(2)P(3) 



(18) 



where we have set up the tetrad frame such that (pi, p 2 , P3) are along (e r , eg, e^) . 

With these expressions we can determine the dispersion relation (15) 
for the Schwarzschiled geometry which can be compared with experimental 
observations. 

3 Polarization dependent time delay 

A radar signal sent across the solar system past the Sun to a planet or 
a spacecraft suffers an additional Shapiro time delay in general relativity 
which has been measured [2] and confirms Einsteins theory to about one 
percent accuracy. In presence of the curvature coupling terms there is an 
additional time delay which we calculate in the following. We consider photon 
trajectories tangential to the gravitating body and in the tetrad frame = 
(p ,Pi, 0,p 3 ). The equation of motion for the (E 1: E 2 ) components is of the 
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p 2 = ±l2\—p 2 3 (20) 



form 

/p 2 + Cl (2(^ 0) -p2 i))+:p 2 3)) \f E A = n 

{ p 2_ Cl{p 2 Q) _ p 2 i)+2p 2 3)) ) [ E J 

(19) 

where C\ = A\(GM/r 3 ). Setting the determinant equal to zero yields the 
dispersion relations for the two propagating modes, 

GM 

We obtain the dispersion relation in terms of the coordinate frame 4-momentum 

by making the transformation p( a ) = e v \ a )k u . The dispersion relations in 
terms of k^ is 

(l-2GM/r)~ 1 k 2 t - (1 - 2GM/r)k 2 - (l/r 2 )k 2 

T 12A^f)=0 (21) 

At the point of inflection of the photon trajectory k r = and one can eval- 
uate the constant angular momentum p^ in terms of the k and the impact 
parameter r as 

fc, = r k t (l - 2GM/rY 1/ \l T 12X™) 1 / 2 (22) 
Substituting (22) in the dispersion relation we get 

(1 -2GM/r)- 1 k 2 - (l-2GM/r)k 2 r -kf\{l-2GM/r)- 1 

, GM , 

We can now solve the dispersion relations for kt in terms of k r to obtain, 
k = k r (l-2GM/r) 

1 (, rl ((l-2GM/r) (l \ (GM_\\\/2 V ' 

\ L r 2 (l-2GM/r )\ L A±[ - rO 3 >> 

where we have defined A± = ±12A. The four-momenta are related to the 
coordinates as k l = dt/ds and k r = dr/ds which enables us to obtain the 
photon trajectory 

dt g lt k t 



=F 12AA; f 2 ^ — (1 - 2GM/r)" 1 = (23) 



dr g rr k r 
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(25) 



Using the relations (25) and (24) we obtain the relation between coordinate 
time t and radius r given by 



t -t = r dr{l ~ 2GM / r ^ i] (2Q ) 

\ L r 2 (l-2GM/r )V 1 A ±y r0 s))' 

Evaluating this integral, we obtain for the time delay 

1 + \/ r /~ r °\ 



At if = r if + 2GMLog{ T 



rf + \jrf-rl 

+ A± ^ ( (!i^ + a^! ) (27) 

The last term is polarization dependent and therefore if one sends a pulse 
of unpolarized radio signal, the two orthogonal polarizations have different 
arrival times at the earth and we should observe a splitting of the pulse in 
time with the time lag between the two polarizations given by 



U-t- = (A + -A_)gy ( "-?- r ° ) + 

ro r f n 

= X 2A^^(\Sin9 f \ + \Sin9 l \) (28) 

where 0, is the angle between the position of the earth and r and Of is 
the angle between the position of the planet /spacecraft and r . The timing 
measurements are accurate to ~ lfisec. When the planet spacecraft is at 
superior conjunction ( SinOi = SinOf = l)and taking M sun = 1.98 x 10 33 gm 
and R sun = 6.95 x 10 5 km, the lack of pulse splitting of the observed radar 
signals to within a jisec accuracy translates to a bound 

A < 1.1 x 10 20 cm 2 . (29) 

This bound is about three orders of magnitude more stringent than the 
one obtained in [8] from the analysis of bending of light by the gravitational 
field of the sun in non-minimally coupled electro-magnetism. 

Since the time is proportional to GM/r 2 this effect is largest for compact 
stars where (M ~ M & and R ~ 10km). Observation of the time delay in 
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binary pulsar PSR B1534+12 [3] yields a value of 6.3±1.3/xsec for the Shapiro 
time delay by the binary companion which has a mass of 1.33M© and radius 
of 10km. The theoretical estimate [4] of this time delay assuming general 
relativity to be the correct theory is 6.6/isec. Taking an upper limit of 1/isec 
as the maximum contribution of the polarization dependent correction to the 
GR time delay we have an get an upper bound on A , 

A < 0.6 x 10 n cm 2 (30) 

This bound is far more stringent than those from the solar ranging exper- 
iments. Note that the constraint on A from pulsars is obtained by taking 
the difference of the travel times of the radio signals with orthogonal po- 
larizations. There are systematic uncertainties in the measurement of total 
travel time from pulsars as discussed in [5]. However these systematic ef- 
fects like the emission time, pulsar velocity, earths velocity and position etc 
are independent of the polarization of the radio signal so these uncertainties 
drop out when we take the difference of travel time between two different 
polarizations. 

For radial photons the four momentum in the inertial frame is p( jU ) = 
(po,Pi, 0, 0) and only non-zero components of the coordinate frame four mo- 
mentum are = (k t ,k r , 0,0). The determinant condition on radial propa- 
gation gives p 2 = and the non-minimal coupling terms have no effect on 
radially propagating electromagnetic waves. 

4 Conclusions 

We find that the curvature coupling of electromagnetic fields give rise to dis- 
persion free photon propagation, where the photon phase and group velocities 
are dependent on the polarization. The best constraints on such couplings 
can be obtained from the timing of electromagnetic waves in the vicinity of 
compact objects. The curvature coupling situation is different from the types 
of couplings where the effect accumulates over the distance of propagation 
and stringent bounds on extra equivalence principle violating couplings are 
put from observations of radio galaxies [7] or gamma-ray bursts [6]. 

In this paper we have taken the view that the classical lagrangian can pos- 
sibly have EEP violating but general covariant terms with large coefficients 
which we bound from observations. If these curvature couplings arise as only 
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as loop corrections to Einteins general relativity [9] and if A ~ M p 2 then 
there are so far no realistic means of observing such couplings in experiments. 
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